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We study the superconducting instability in systems with long but finite ranged, attractive, 
pairing interactions. We show that such long-ranged superconductors exhibit a new class of fluc- 
tuations in which the internal structure of the Cooper pair wave function is soft, and thus lead 
to "pseudogap" behavior in which the actual transition temperature is greatly depressed from its 
mean field value. These fluctuations are not phase fluctuations of the standard superconducting 
order parameter, and lead to a highly unusual Ginzburg-Landau description. We suggest that the 
crossover between the BCS limit of a short-ranged attraction and our problem is of interest in the 
context of superconductivity in the underdoped cuprates. 
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I. INTRODUCTION 

Superconductivity is a subject that has attracted prolonged interest in the condensed matter comniunity. Much of 
our current understanding of this subject is based on the standard Bardeen-Cooper-Schrieffer (BCS)tl theory, which 
has enjoyed great success when applied to conventional superconductors. The discovery of the "high-Tc" cuprate 
superconductors opened a new chapter in the study of superconductivity. It is now well understood that just like 
in the BCS theory, electrons in cuprate superconductors form Cooper pairs. Moreover, it has been established that 
the wave functions of these Cooper pairs have predominantly d-wave symmetry. However, the origin of the force 
that leads to Cooper pairing in these systems is not yet understood and, experimentally, significant deviations from 
the BCS framework (which we take tO|-jnclude a Fermi liquid description of the, normal state) have been found near 
and above the transition temperaturep especially in the underdoped cuprates.B Theoretically, there is a continuum 
of proposals for understanding these anomalies which tange from conceptually modest (if_practically far reaching) 
enhancements of the BCS framework (spin fliictuations,Q the crossover to pre-formed pairs,Q etc), to the influenajpOj 
various hypothesized quantum critical pointsEl on to spin charge separation and non- Fermi liquid normal states.au' I 
Finally there are scenarios based on self-organized dimensional reduction (stripes) which are very different in spirit.E 

This paper, also motivated by the physics of the cuprates, is firmly in the conceptually modest camp. We examine a 
new crossover out of the BCS corner, this time to a superconductor with long range pairing interactions, i.e. in which 
the range of the pairing potential L is much longer than the coherence length (L 3> deduced from measurements 
of the (anisotropic) gap in the quasiparticle spectrum. We remind the reader that in the BCS theory, the phonon- 
mediated, retarded attraction between electrons is often modeled as a short-range instantaneous attractive potential. 
As a matter of fact, it is often modeled as a consta nt . at tractive scattering potential in momentum (k) space, which 
corresponds to a (5-function attraction in real space.t2rtil This is an excellent approximation because in conventional 
superconductors the coherence length ^ is of order ~ lOOOA, while the phonon-mediated attraction is a local (but 
retarded) on-site interaction with range of order the lattice spacing, and therefore the range is effectively zero compared 
to the size of the Cooper pairs. In such a short range model, the elementary excitations are quasiparticles with a 
(k independent) excitation gap A; the only collective mode below 2A is the linear Goldstone mode, whose energy is 
pushed up to the plasmon energy in real systems due to the epdstence of the long-range Coulomb interaction. The 
classic BCS limit also involves the assumption of weak couplingi3 whereupon thermally excited quasiparticles control 
the largely mean-field superconducting transition temperature T^, giving rise to a universal ratio 2A/Tc « 3.5. 

Our interest in what happens when the range of the interaction is no longer the lattice spacing but instead crosses 
the (suitably defined, see below) coherence length arises from four considerations. 

First, ^ is much shorter in the cuprates than in conventional superconductors; in fact ^ is of order five lattice 
spacings in the a-b plane, and less than one lattice spacing along the c-axis. Therefore even if it turns out that the 
effective attractive interaction that gives rise to pairing is of the order of a lattice spacing, L and ^ will be comparable. 
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Second, the physics of the pseudogap regime is most simply explained by invoking a suppression of by fluctuations 
that are small in the BCS limit. Quite generally the order parameter (or pair wavefunction) involves a relative piece 
and a center of mass piece, and the latter is what enters the standard Ginzburg-Landau theory. The simplest possibility 
is that the phase of the center of mass piece fluctuates strongly and this has been suggested in the context of the 
cuprates from two different standpoints, asJic physics of strongly coupled pre- formed pairsO or that of a low superfluid 
density coming from a doped Mott state.EJ A more novel possibility is that the relative pair wavefunction has soft 
fluctuations, and as we demonstrate in this paper, these arise in the limit of a long range attraction. 

Third, it appears that the model we study here is relevant to some of the proposed mechanisms for cuprate 
supercoiidjjjctivity. The first example is the inter-layer pair hopping mechanism proposed by Anderson and co- 
workers.Bo In this theory, pairing is induced by coherent pair hopping along the c-axis, pipvided that single electron 
hopping is frustrated by the non- Fermi liquid nature of the normal state. It was emphasizecuJ that high comes from 
the fact that the matrix element for pair hopping is diagonal in k space, although the symmetry of the order parameter 
is determined by the in-plane pairing potential. Mathematically, the pair hopping term is equivalent to an attractive 
potential off-diagonal in layer index, even though its physical origin is kinetic energy along the c-axis; in particular, 
being diagonal (or a 5-function) in k space, it corresponds to a pairing potential that has an infinite range in real 
space. A second example of this is the spin-fluctuation mechanism proposed by Pines, Scalapino, and co-workers.a 
In this theory, the range of the effective pairing interaction mediated by antiferromagnetic spin fluctuations is the 
spin-spin correlation length; as one reduces the doping level and approaches the antiferromagnetic instability from 
the overdoped side (where the system is more BCS-like), the spin-spin correlation length and therefore the range of 
the pairing interaction increases; it can become very large in the underdoped region, and diverges when approaching 
the antiferromagnetic instability. 

Finally, given the importance of the subject of superconductivity, we feel the situation we study here (which is 
opposite to the familiar short-range attraction limit) is interesting in its own right. Eiiep if it proves not to be of 
particular relevance to superconductivity in cuprates, it may become relevant elsewhere.lij 

Our main results may be summarized as following. We find in addition to the usual quasiparticle excitations, the 
system supports collective modes whose energies are significantly lower than the quasiparticle gap A, when L ~ ^; in 
the limit that i S> ^, the lowest energy collective mode gap becomes much lower th&n A, and the number of collective 
modes below 2A becomes large. At finite temperatures, due to the thermal fluctuations of these collective modes, 
the transition temperature Tc becomes significantly lower than the mean field transition temperature, T^^^ (which is 
controlled by the thermally excited quasiparticles and does not know about the collective modes). For < T < T*^^, 
the electrons are paired (in momentum space) but there is no long-range phase coherence; thus the system exhibits 
pseudogap behavior. Also the effective Ginzburg-Landau theory needs to be modified in this case; in addition to the 
center of mass degree of freedom of the superconducting order parameter (or Cooper pairs), which is the only degree 
of freedom taken into account in the original Ginsburg-Landau description, we also need to take into account the 
fiuctuations of the relative or internal degrees of freedom of the order parameter. As a consequence spatial gradients 
couple to the internal pair structure in non-trivial ways, and although we have not done the computation yet, we 
suspect that the critical field Hc2 will be governed by a coherence length that is distinct from the one derived from 
the gap. 

In the main part of this paper, we will take the attitude that this is an interesting model system to study, and 
work out various properties of such a model, especially those that are different from the short-range limit of the BCS 
theory. Discussion of the possible relevance of our considerations to the cuprates will be reserved to the Summary 
section. 

The rest of the paper is organized as following. In Section || we revisit the Cooper problem, and keep the finite 
range of the attractive potential explicit. We note some problems involved in taking the limit of long-ranged,potentials 
and specify the precise nature of the limit we consider in this paper. We find that in the limit L ^ ^Jlj there are 
a large number of bound state solutions, whose binding energies are comparable to the ground state. These bound 
states correspond loosely to the collective modes in the many-body problem. In section |in| we study the BCS reduced 
Hamiltonian with a finite range attractive force, solving for both the mean-field ground state, and the zero momentum 
collective modes using the linearized equations of motion for the order parameter. In section ^ we develop a Ginzburg- 
Landau effective theory from our model, using the functional integral method. We find that it is necessary to keep 
track of the fluctuations of both the internal and center of mass degrees of freedom of the superconducting order 
parameter; in particular, in the limit i 3> ^, Tc is controlled by the fluctuations of the internal degrees of the order 
parameter, and is much lower than T*^^. In Section ^ we calculate the reduction of from T^^^ due to the thermal 
fluctuations of the collective modes, using the Ginzburg-Landau theory developed earlier. In Section ^ we discuss 
the possible relevance of our work to the cuprate superconductors and mention some connections to other work in the 
literature. 
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II. THE COOPER PROBLEM 



We begin by studying the Cooper problem, in which the wave function of a Cooper pair takes the form: 

|vl/>= J2 gik)clcl^^\^o), (2.1) 

0<ek<Ea 

where \'^o) is the fiUed Fermi sea, — vpik — kp) is the single electron energy measured from the Fermi level, and 
Ec <C Ep (Fermi energy) is an energy cutoff. In phonon mediated attraction models, Ec is usually taken to be the 
Debye energy. In this paper we treat Ec as a parameter that can be varied in our effective model. The Hamiltonian (for 
both the Cooper problem considered here and the many-body problem in the BCS reduced Hamiltonian approximation 
considered in the following section) takes the form 

-f^ = X] "^fe^aCka - X! ^k-k'|4|clkxC-k'iCk'T. (2.2) 
k,(T kk' 

Here 

Vq^jJ d2rF(r)e-'i'" (2.3) 
is the Fourier transform of an attractiv e po tential V{r); A is the area of the system.lll It is understood that the 



summation of the second term of Eq. (2^) is restricted to states with |efc|,|efc'| < Ec. For simphcity, we assume 
Vq and V{r) are positive definite, so the ground state for the Cooper problem has s-wave symmetry, and the many- 
electron system is an s-wave superconductor. We will comment on the generalization to a d-wave superconductor in 
the concluding section. 

The Schroedinger equation that g(k) satisfies is 

Eg{k) = 2efc5(k) - ^ %_k'|5(k'). (2.4) 



where E is the eigenenergy. For a bound state, we must have E < 0. For a general Vq, the integral equation (2.4) is 
difficult to solve. Cooper solved the short-range limit of this problem. In this limit, Vq ~ V is a constant in momentum 
space, and therefore a 5-function in real space. It was found that there exists only one bound state solution, with 
s-wave symmetijy; for weak coupling {N{0)V ^ 1, A^(0) being the density of states at the Fermi energy), the energy 
of this state isK 

E ^ ~2Ece-^/^^'>^^ . (2.5) 

This well known dependence of the binding energy on V is extremely singular. For finite but very short-range attractive 
potentials, there may be more bound state solutions (probably in other angular momentum channels). However, due 
to the singular dependence of the binding energy on the pairing potential, the binding energies of these states will be 
much smaller than the ground state, even if the effective pairing potential only changes slowly from one channel to 
another. 

Here we consider the opposite limit, in which V{r) has a long but finite range, L. Although the explicit form of 
V{r) is unimportant to our basic conclusions, for concreteness we assume it has a Gaussian form: 

F(^) = Foe-'-'/2^', (2.6) 

and therefore 

Vq = ^Foe-^^^/^. (2.7) 

We are interested in following the evolution of the system as L is increased and in the regime where L is large, 
L > ^. One possibility is to consider potentials of the "Kac type" familiar from statistical mechanics, where the 
range is increased while keeping the integrated strength of the potential fixed, thereby achieving an unproblematic 
thermodynanaic limit. A second possibility, perhaps more appropriate to attractive interactions generated by the 
system itself,t3 is to keep the interaction of fixed magnitude (Vq) but to restrict its operation to an increasingly 
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narrow shell around the Fermi surface. We will choose the second course here. An easy estimate shows that we need 
to pick the cutoff energy 



-D 



(2.8) 



in order to keep the potential energy per particle finite. This also has the advantage that, by construction, is suppresSjes, 
the tendency to phase separation that would otherwise be a complication with attractive long range interactions .tS 
There is however, a final caveat. For technical reasons, starting in the next paragraph, we will often employ a gradient 
expansion in momentum space. This is not quite consistent with the cutoff procedure but for the results of interest 
we will find that the precise value of the cutoff will enter weakly, leading us to believe that a better set of calculations 
will not alter our general scenario. 

Since V|k-k'| goes to zero rapidly for |k — k'| > 1/L, if g(k) var ies slowly on the scale of 1/L, we may perform a 
gradient expansion for g(\i) in k-space, in the last term of Eq. ( |2.4D : 



E^k-k'|3(k')=EMk-k'| 



5(k)+Vk5(k)-(k'-k) + i^ 



2 ^ dkadK ^ 



(E + J Ik - k?%-k'| j vLg(k) 



^05 (k) 



A. 

2L2 



V|5(k), 



in which we have neglected higher gradient terms. Thus within the gradient expansion, Eq. (2.4) reduces to 

1 



2M 



V^ff(k) + (2efe - yo)ff(k) = i?ff(k). 



(2.9) 



(2.10) 



where M — I? jV^. This differential equation (2.10) is identical to the Schroedinger equation of a particle confined in 
an annulus (with some unknown, but calculable, boundary condition) with inner and outer radii kp and kp -\- Ec/vp 
respe ctivel y, experiencing a "potential" U{k) = — Vb + 2efe = — Vb + 2vp{k — kp). It is interesting to note that in 
Eq. ( |2.10| ), the "kinetic energy" term comes from the two-body interaction in the original Hamiltonian, while the 
"potential" term actually is the original kinetic energy term. 



For large L, the "mass" M is large, and we may use WKB approximation to analyze Eq. (2.10). The following 
conclusions follow straightforwardly: 

i) The ground state is in the s-wave channel, with energy 



Eq 



-Vo + c{ 



Vp\ 



L 



(2.11) 



Here c is a constant of order 1; it is approximately (^^)^^'^ for hard wall boundary conditions. We find the binding 
energy is essentially the depth of the attractive two-body potential Vb ; this dependence is much less singular than the 
short-range limit defined above. Also the cutoff Ec does not enter explicitly.E3 At this point we define a "coherence 
length" 

e = 2vpI^\Eo\ » 2vp/^Vo, (2.12) 

in analogy to the BCS coherence length at zero temperature: Co = "v p / {tt /S.{Q)) , where A(0) is the quasiparticle 
gap at zero temperature. What we have in mind is that the binding energy in the Cooper problem will be the same 
as twice the quasiparticle gap in the BCS problem. This is not the case in short-range models; however, as we will 
see later, it is indeed true in the present case. We need to emphasize however, that <^ is not the size of the Cooper 
pair wave function t in the two-body problem studied in this section (this is the case even in short-range models); the 
latter can be esti mate d easily from Eq. (2.10). The size of the ground state wave function in momentum space for 
the Hamiltonian ( 2.10| ) is 



Ak 



1 



vpM 



1/3 



KvpL-^' 



1/3 



(2.13) 



thus I ^ 1/Ak ^ (vpL^ /Voy^'^; it actually increases with L; however the sublinear dependence means for large L 
we do have £ L. This points to some ambiguities in what is meant by the coherence length in our problem and 
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presumably spatial gradients may be governed by a different coherence length than the gap coherence length. We 
note that this possibility has been raised on completely different grounds in Ref. |^. 

In order for the gradient expansion performed above to be valid, we must have Ak > 1/L, i.e., (/(k) of the ground 
state varying slowly over the scale 1/L, which is the range of Vq. This leads to the condition 

L > 2vf/7:Vo « (2.14) 

as advertised earlier. It turns out that throughout this paper the gradient expansion in momentum space is the key 
technique that allows explicit calculations of various physical quantities to be made, and the above condition defines 
the range of its validity. For a given potential with fixed Vq and L, on e may also interpolate between the short-range 
and long-range limits by varying vp- In the following we assume Eq. ( 2.14| ) is satisfied, and pay particular attention 
to the limit L 3> 

ii) The number of bound states in the s-wave (/ = 0) channel is A^'o ~ ^^^^ • 

iii) The largest angular momentum channel that supports a bound state has angular momentum Imax ~ V2kpL. 

iv) The total number of bound states is Ntot = ^i^i ^ f ^f^^/C- 

v) The energy spacing between the ground state and the first excited state is AE w l-^^l' much 
lower than the binding energy of the ground state itself, in contrast to the short range case. 

The gradient expansion does not apply to all bound states, even if L ^ ^. One can show that it only applies to 
states with energy measured from the ground state (or bottom of the potential well) AE < VbC^/i^. Nevertheless 
the estimate of number of bound states in various channels should be qualitatively correct. 

These large number of low-energy bound states correspond loosely to the low-energy collective modes in the many 
body problem studied in the following sections, and they lead to a significant reduction of from its mean-field value 
(whose scale is set by the gap which is the same as the Cooper pair binding energy here, T*^^ ~ A ~ Vq). This 
may be understood heuristically in the following way. When the temperature T is reduced to T^^^ , electrons start 
to form Cooper pairs. In the short-range BCS model, there is only one way to form a Cooper pair, thus all pairs 
condense into one state immediately and long-range coherence is established. In the model studied here however, 
there are many different ways to form Cooper pairs and their binding energies are comparable; thus the system needs 
to go to much lower temperature for all the Cooper pairs to condense into the ground state, hence Tc <C T*^^. 
For Tc < T < T^'^, the electrons are paired (and therefore gapped), but there is no long-range superconducting 
coherence, and the system exhibit pseudogap behavior. 



III. MEAN FIELD THEORY AND COLLECTIVE MODES IN THE BCS REDUCED HAMILTONIAN 



In this section we study the BCS reduced Hamiltonian, Eq. (2.2). Instead of solvi ng t he two-body problem in 



the above section, here we study the many-body problem. The difference between Eq. (2.2) and the full many-body 



problem is that in Eq. (2.2), only pairs of electrons with opposite spin and total momentum zero scatter each other, 
(we return to the full problem, i.e. when the sum over the interacting momenta is constrained only by momentum 
conservation in the next section.) This leads to a very special property, that an unpaired electron never gets scattered, 
and thus the number of unpaired electron (0 or 1) for any k and spin orientation is a conserved quantity. The ground 
state for this Hamiltonian is in the subspace in which no such unpaired electron exists, namely the pair of single 
electron states k | and — k J, arc cither both occupied or both empty. In this subspace, one may map the problem 
onto a quantum pseudospin problem.cj We introduce a pseudospin-1/2 operator for each k: 

= (cLt^kT + clkiC-ki - l)/2, 

■^k 

'5k=c-kiCkT. (3.1) 



^^ckT -r L,_ki' 
^k — '-kT''-ki' 



In this mapping, the pseudospin is pointing up when a pair of single electron states are occupied, and down when 
they are empty; pair creation/annihilation operators map onto spin-flip operators. The Hamiltonian now takes the 
form (up to a constant): 

H^J2 2ek^k - ^ E Mk-k'|(^,[5k, + S+S^) = J2 2ek5£ - J2 %-k'l iS^S^> + Sl,Sl), (3.2) 

k kk' k kk' 

i.e., the kinetic energy maps onto a Zeeman field that depends linearly on (and changes sign at the Fermi surface), 
which couples to the z component of the pseudospins; the pairing interaction maps onto a ferromagnetic coupling 
among the xy components of the pseudospins. 
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In the mean field approximation, one replaces the ferromagnetic coupling among the xy components of the spins 
by an average field: 



-ffji/F = — ^ • Sk, (3.3) 

k 

where satisfies the self-consistency equation: 

=-2eki + 2^1^k-k'|(S;^'). (3.4) 

k' 

Here stands for the xy components of Sk. In the ground state, Sk points in the direction of B^; the elementary 
excitations are single pseudospin flips, which corresponds to a pair of Bogliubov quasiparticles with opposite momenta 
and (real) spin.Ej 

In the short-range limit L —^ 0, the range of V|k_k'|i ^/L, diverges in k space. Thus all the pseudospins are equally 
coupled, no matter how far away they are in k space. In this limit, the mean field approximation becomes exact; 
the only collective mode is a zero mode corresponds to-4he global rotation of all the pseudospins along the z-axis, 
refiecting the broken XY symmetry of the ground state.cil All other excitations in this subspace may be described by 
pseudospin flips, or pairs of quasiparticles. 

The situation, however, becomes quite different, when L becomes large. In this case V|k-k'| becomes short-ranged 
in k space; one may actually divide the k space into blocks of size x ■^; within each block all the spins are strongly 
coupled and are effectively locked into a very big single spin; for different blocks, however, couplings are restricted 
to neighboring blocks. In such an XY ferromagnct with "short range" couplings among the blocked spins, there 
are (pseudo)spin-wave like collective excitations in addition to single pseudospin flips, whose energies can become 
significantly lower than pseudospin flips for large L, as the number of blocks increases as ^ L^. 



We begin by analyzing the mean field equation, Eq. (3^), at zero temperature. Without losing generality, we may 
assume B|^ is in the x — z plane: 

B|^ = -2ekz + 2AkX, (3.5) 

and the self-consistency equation for Ak is 

k' ^ k' V ^k' + ^k' 



which is the familiar BCS gap equation.E2[Eil In the limit L 0, Vik-k'i = ^ becomes independent of |k— k'|, so does 



Ak = A, and Eq. (3.6) reduces to 



In the weak coupling l imit , its solution is A = 2i?ce~^/^*^°^^. This gap is much larger than the binding energy in the 
Cooper problem, Eq. (2.5), due to the factor of two difference in the exponential. 

The situation becomes very different in the opposite limit, that L becomes large. In this case V|k_k'| becomes 
short-ranged in k space; thus on the right hand side of Eq. ( [3.6D , we may replace Ak' and Ck' by Ak and ek to 
first approximation, provided that they vary slowly on the scale of 1/ L. In this approximation one finds in the limit 
L — > oo, 



Ak = y^-4. (3.8) 

for |ek| < Vb/2, and Ak = otherwise. For large but finite L, Ak is nonzero but exponentially small for |ek| > Vb/2. 
Thus Ak is k dependent^ and reaches its maximum for fc = fc^?, in which case 

Afc, = (3.9) 

We find the quasiparticle gap and the binding energy are indeed set by the same energy scale in the Cooper problem, 
as advertised earlier; in particular, the gap for creating a pair of quasiparticles on the Fermi surface, 2Afey — Vb, 
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is exactly the Cooper pair binding energy, in the Hmit L — s- oo. Following standard convention jl^ we introduce the 
coherence length 



(3.10) 



This definition indeed matches that of ^ in the Cooper problem Eq. 2.12| , in the limit that L is large. And one can 
easily show that the large L approximation is valid for 

L>C (3.11) 
Within the mean field theory, the elementary excitations are the Bogliubov quasiparticles, with the standard spectrum 



Eh — 



At finite temperature, the self-consistent mean field equation becomes 



tanh 



(3.12) 



(3.13) 



where /3 = is the inverse temperature (we set the Boltzmann constant ks to be 1). We can use the above equation 



to de term ine the mean-field transition temperature T*^^. 
Eq. (3.13) reduces to 



In the L — > oo limit, it is particularly simple; in this case 



1 



Vo ^ , /3i?k 

tanh . 

2Ek 2 



(3.14) 



Using the facts that as T ^ T^^ , A 
the strongest, we find 



k — !■ and i?k eic) a-nd focusing on |k| ~ kp where the pairing instability is 

(3.15) 

in this limit. Thus the ratio 2Af^p /T^^^ — 4, which is slightly bigger than the BCS ratio of 3.5. 

There are a couple of new features of the mean-field solution in the large L regime, that deserve some further 
discussion, (i) Unlike the case of short-range attraction originally considered by BCS, where the gap A is a constant 
in momentum space at all temperatures (within the cutoff where the pairing interaction is nonzero), here Ak has 
strong momentum dependence, being m axim um at the Fermi surface and decreasing as one moves away from it. At 
T = 0, the momentum dependence Eq. (3.8) is such that the quasiparticle dispersion is essentially flat near the Fermi 
surface, for sufficiently large L. (ii) The temperature T not only affects th e ove rall scale of the gap, but also its 
momentum dependence. Taken at face value, one would conclude from Eq. (3.14) that T*^^ depends on ek in the 
following way: 



MF 



2ek/Fo 



tanh"'(2ek/Fo) 



< T, 



MF 



(3.16) 



What this really means, of course, is that the gap Ak remains exponentially small for ^^^^(ek) ^ T < T^^^ . Thus 
in some sense T*^^ is a "local" property in the momentum space, increasing monotonically, and scaling roughly with 
the size of local gap Ak(T — 0) unless Ak(r = 0) is very small. 

In the following we go beyond mean field theory and study the collective excitations of the system at T = 0. To do 
that, we study the equations of motion for Sk": 



where 



^ = -i[Sk,H]=SkxBk, 
at 



Bk = -2ek5 + 2^14k-k'|S^,. 



(3.17) 



(3.18) 
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Write Bk = + 5Bk and Sk = + (5Sk, where is the expectation value of Sk in the mean field ground state, 
and linearizing the equation of motion by neglecting terms proportional to (5Sk x (5Bk, we obtain 

^ = 5Sk X BO + 8° X 5Bk. (3.19) 

(5Sk should be perpendicular to Sk- If Sk is in the x — z plane, we may assume 

5Sk = 5Sly + 5^11 ek, (3.20) 
where Ck is in the x ^ z plane, but perpendicular to S^. Thus Eq. ( 3.19D reduces to 

-^ = -5Blcose^-Bi6Sl (3.21) 

^ = Bi5Sl \5Bl (3.22) 
Here is the angle between Bj^ and the z direction, and 

5BI = 2^%_k'|'55j|,cos0k', (3.23) 

k' 

<5i3^ = 2^14k-k'|'55^'. (3.24) 

k' 

Thus for a mode with frequency lj and 6S^ oc 0k, it must satisfy 

- cc;20k = -{Blf<t>^. + Mk-k-l-Z-k' + cos 0k ^k-k'|S^' cos(?k'0k' 

k' k' 

- cos 61k ^ V|k_k'|Mk'-k"| cos6'k'0k". (3.25) 

k'k" 

Performing a gradient expansion similar in spirit to the previous section, we obtain 

c^'0k = f/k0k - \Bi ^ |k - k'|2^k-k'|V20k - ^(J] |k - k'|2^k-k'|)V2(cOS0kS£(/'k) 
k' k' 

+ cos^ki^ V|k-k'|V|k'-k"| cos6'k'(k" - k)] • Vk<?!>k 

k'k" 

+ icos0k[^%-k'|%-k"|Cos0k'(fcM" -fc-)]5M5-<^k + ---, (3.26) 

^ k'k" 

where the ellipses stand for terms that involve higher gradients which we neglect. Here 

J/k = {Blf - Si:^ ^^k-k'l - cos^ 0k(i325] V^k-k'i - E %-k'|^k'-k"|) - COS0, 

k' k' k'k" 

= {Bl - cos" e,,Vo){Bl ~ K,) - cosflfc. (3.27) 

In the limit that L becomes very large, one finds C/k = for |ek| < Vo, and f/k = (fk — Vb)^ otherwise. 

For low energy modes, we expect 0k to be centered near k = kp. Thus in a somewhat crude approximation, we 
may neglect the gradient terms in Eq. (3.26) proportional to cos6'k, as cos 0k — > as A: ^ kp] and assume C/k to 
become infinite for |ek| > Vq. Within this approximation (which should give qualitatively correct results), we obtain 
the collective mode frequencies to be 
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where m and n are "momentum" -like quantum numbers along and perpendicular to the Fermi surface respectively. 
So the mode frequency is linear in these "momenta" , as one would expect for the spin wave spectrum in an XY 
ferromagnet. It is clear that the energies of these collective modes become much lower than those of single pseudospin 
flips (corresponding to quasiparticle excitations), when L is large. Since our discussion in this section is restricted to 
the BCS reduced Hamiltonian in which only pairs with zero total momentum interact with each other, these modes 
are exciton-like collective modes at zero momentum; they will acquire the usual dispersion of exciton modes at finite 
momentum. The term "exciton" indicates that the zero momentum excitation consists of a pair of quasiparticles 
bound into a pair state different from that of the condensate which is why the state lies in the gap between the ground 
state and the two quasiparticle continuum. The availability of such states (as we already saw in the previous section) 
is due to the long range of the interaction. 



IV. GINZBURG-LANDAU THEORY AND ELECTROMAGNETIC RESPONSE 



In this section we use the functional integral formalism to derive the effective Ginzburg-Landau theory near their 
jiMF Jqj. superconductors with finite range attractive interactions, and use it to derive their analog of the London 
equation. 

Let us consider the Hamiltonian 

H^f + V^ ^(efc - M)cLck. - J dxdx'F(|x- x'|)*j(x)*j(x')*i(x')*TW- (4.1) 

k(7 

Here V{x) > represents an attractive interaction. We will primariliL-be interested in singlet pairing in the ground 
state; so we neglect interactions between electrons with the same spin.L3 While the interaction is written in a pairwise 
form, it is understood that a cutoff exists in momentum space so that only electrons that are close enough to the Fermi 
surface interact with each other, in the same manner as in previous sections. The difference here is that the interaction 
is no longer restricted to pairs of electrons with zero total momentum; thus the above Hamiltonian represents the full 
many-body problem, albeit with the parallel spin interaction ignored. 

One may also describe the system using an Euclidean action in terms of Grassman variables: 

S[i','^] = So[i','^]- dr J (ixdx'y(|x - x'|)¥|(x, t)¥x(x', t)*! (x, t)*| (x, r), (4.2) 
where So is the action for free electrons, and r is the imaginary time. The partition function is 

i?*Z?^e-^[*'*l. (4.3) 



We now decouple the quartic term in S by introducing a pair of Hubbard-Stratonovich fields A(x, x', r) and A(x, x', r), 
which will become the superconducting order parameter in the sequel: 

5[*,*,A,A] = So- dr j dx<ix'[A(x,x',r)l'|(x,T)*|(x',T) + A(x,x',r)«'|(x',r)^'|(x,T)] 



J 

With this decoupling, the fermionic action becomes quadratic, and can be integrated out, after which we obtain an 
effective action in terms of the order parameter A(x, x', r): 

5e[A,A]=/ drdxdx' ' ' ' , ^logZ[A,A], (4.5) 

Jo K (|X — X |j 

where 

Z[A ^] = y D^D^e~^°^'^^^^^o ■'"'^'^^^^''''''^ (4 6) 

The mean- field solution corresponds to the saddle point of 5*6 [A, A]: 
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SSe[A,A] 

(5A(x,x') 



A=A, = 



A,(x,x') 5 log Z[ A, A] 



V^(x-x') 5A(x,x') 



A=A, 



|^-(*,(x'm(x)>A.^O, 



(4.7) 



where ()as stands for quantum and thermal averaging in the presence of the pairing field As(x,x'). Here we have 
assumed a static saddle point so that As has no r dependence. F urth er assuming that at t he s addle point, As(x, x') — 
As(x— x') is translationally invariant, it is easy to show that Eq. (4.7) is equivalent to Eq. (3.6), upon the identification 



A 



dxe 



As(x) 



(4.8) 



The functional integral formalism can be used to derive the effective Giiizburg-Landau free energy, in the vicinity 
of . This has been done for the short-range attractive interactions.Ej Here we use it to derive the appropriate 
Ginzburg-Landau free energy for finite range attractive interactions. 

Our starting point is the effective action, Eq. ( [4.51 ). Near T*^^, we may make two simplifications: i) We may 
neglect the t dependence of A as we expect the thermal fluctuations to dominate the quantum fluctuations; ii) We 
may expand Se in powers of A. The quadratic terms take the form 



5(2)[A,A]=/3 j dxdyi^^^^y dxidyidx2dy20(xi,yi;x2,y2)A(xi,yi)A(x2,y2), 



(4.9) 



where 



(xi,yi;x2,y2) 



52 logZ 



^A(xi,yi)(5A(x2,y2) 



|A=0 



dTi / dT2(^'x(y2,r2)1't(x2,r2)*t(xi,Ti)^'i(yi,Ti))c 

Jo 

= dn dr2Go(y2 - yi;T2 - ri)Go(x2 - xi;t2 - Ti) 
Jo Jo 

= ^Go(y2 -yi;ia;„)Go(x2 -xi;-ia;„). 



(4.10) 



Here Go is the normal state (non- interacting) single electron Green's function, uj^s are fermion Masubara frequencies, 
and ()c stands for connected contractions in the average. 

We now introduce "center of mass" and "relative" coordinates for the order parameter A(x, y): 



R=(x + y)/2, r = y-x, 

and Fourier transform with respect to the relative coordinate r: 



A(R, k) = y" dre-^''"'A(R - r/2, R + r/2). 



In a uniform (R independent) conflguration, A(R , k) becomes Ak- 
In terms of A(R, k), the first term of Eq. (4.9) becomes 



(3 / dxdy 



|A(x,y)p 
nix-y|) 



(3 dR dr 



|A(R-f,R+§)p 
V{r) 



A. 



where 



j dRY. ^(^' ki)A(R, k2)F(|ki - k2|), 

•' kik2 

F(|ki -k2|) = j dre'(''^^''^>7y(r). 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



This expression is problematic as it stands, for 1/V{r) does not, in a strict sense, possess a Fourier transform in a 
truly infinite system. However, it can be defined in a large but finite-size system with periodic boundary conditions; 
and as we will see in a few lines, this is only an intermediate expression and can be regulated as such and the final 
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answer will be sensible even as the regulator is removed. It is easy to see that such a regulated F{k) must be peaked 
near fc = 0, and therefore short-ranged in k space. Thus we may perform a gradient expansion for A(R, k) in k space, 
as in previous sections. Introducing k = (ki + k2)/2, k' = k2 — ki, we obtain 

^ fdHYl ^(^^ ki)A(R, k2)F(|ki - k2|) 

= 4 / dRE(l^(I^'k)p~^|VkA(R,k)p + ...)F(fc') 
kk' 

f dR[Si|A(R,k)|2+S2|VkA(R,k)|2 + •••], (4.15) 



A 



where 



'^'-jj:n'')-^y ("6) 

k 



k 

1 

B 



k 



For V{r) ~ Vqc ^^/^^^^ -^e have 



Bi - ^, (4.18) 

Clearly, the coefficient B2 controls the fluctuations of A in the k space, which describes the internal degrees of freedom 
for the pairing amplitude; the larger the interaction range L is, the softer such fluctuations are. On the other hand, 
in the short range limit L 0, such fluctuations get completely suppressed, and the order parameter A depends only 
of the center oLpiass coordinate R, and we recover the standard Ginzburg-Landau theory in terms of A(R) with no 
k dependenccEj It is clear however, we need to keep full k dependence of A in the present case. 

In the first term of Eq. (4.9), there is no mechanism to control the spatial (R) fluctuations of A. Such fluctuations 
are controlled by the second term, which we now turn to. The second term in Eq. (O) takes the form 

- / rfx,dy,dx2dy2Q(x„ y, ; X2, y2) A(x„ y,) A(X2 , y2) 

= -jJ2 [ o«Rio«R2A(Ri, k)A(R2, k)P(R2 - Ri, k), (4.20) 

k 

where 

P(R2-Ri,k)-i Y Go(ki,zc.„)Go(2k + ki,-*L^„)e^''''+'''^^-^''="''^\ (4-21) 



A 

kl jlLJn 



and 



Go(k,zc^„) = . . (4.22) 

lU!n - Ck 

In this term there is coupling between A's with different R's, but no coupling between A's with different k's. Thus 
the first and second quadratic terms in the action control the internal (k) and spatial (R) fluctuations of the order 
parameter A respectively. 

The fluctuations of the overall magnitude of the order parameter are controlled by higher order terms in the effective 
action. As usual we may stop at the quartic term near T^^^ , and neglect terms that involve the spatial gradient of 
the order parameter at this order .0 For the present problem, we obtain (in a way similar to Ref. 24 ) 

Si'' - AE / dR|A(R,k)|^(^^^^i^^^ - 1 ). (4.23) 

SA^J I V , ;i V ^3 2e2cosh2(/3ek/2)^ ^ ' 
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Now we perform a spatial gradient expansion in R, for the second quadratic term, Eq. (4.20). To do that, we first 
go to momentum space, and define 



A(Q,k) = ydRe-^'5-^A(R,k). 
In terms of A(Q,k), the second quadratic term reads 

l^D(Q,k)|A(Q,k)|^ 



A2 



Q,k 



where 



D(Q, k) = ^ Go(Q/2 - k, iu;„)Go(Q/2 + k, -iu,,) 
_ /3tanh(/3eQ/2-k/2) +tanh(/3eQ/2+k/2) 



Expanding for small Q, we obtain 

D(Q,k)«i?(0,k) 



eQ/2-k + eQ/2+k 



13^1 sinh(/3ek/2) 



32ekCOsh-'(/3ek/2) 
Fourier transforming back to real (R) space, we obtain 



(Q-k)= 



i^i?(0,k) /dR|A(R,k)p + ^5] 



A 



1 



P^'vl sinh(/3ek/2) 
^^32ekCOsh3(/3ek/2) 



dR|k- VA(R,k)|^ 



Thus putting all these terms together, and keeping leading gradient terms only, we obtain 
5e[A,A] = ^^ / dR{[Bi-C(/3,k)]|A(R,k)p+B2|VkA(R,k)p 



|k.VA(R,k)|2 + [/(fc)|A(R,k)|n, 



where 



Bi 

B2 



1/^0, 

l/(2V^oi'), 



C(/3,fc) =tanh^/(2efe), 
1 sinh ^ 



2TOk 32efe cosh^ ^ ' 



f/(fc) = 



tanh 



84 



16e2 cosh^ ^ ■ 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 

(4.30) 
(4.31) 

(4.32) 
(4.33) 

(4.34) 



Eq. (4.29) is one of the central results of this paper; it is the basis for the calculation of reduction of Tc from T*^^ 
due to thermal fluctuations of collective modes in the next section. 

Thus far in our discussion we have assumed a uniform space with no background electtomagnetic field. In the 
presence of a weak and slowly varyi ng ba ckground magnetic field, using standard argumentsEJ one finds that the only 
modification one needs to make in (4.29) is to replace V by V + 2ieA(R)/c. The supercurrent is 



2ieA 



)A] 



■}• 



(4.35) 



^From this we may obtain the equivalent of London's equation: 
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J(-) ^ ^ E ■ A(H)1 -i^E ^)A(K,, 

k k 



from which the expression for penetration depth follows: 

(4.37) 

Here d is the spacing between layers. These are, of course, mean-field results; fluctuations have been left out at this 
level. And they apply only near T*^^, even at the mean-field level. 



V. REDUCTION OF Tc DUE TO THERMAL FLUCTUATIONS OF COLLECTIVE MODES AND 

PSEUDOGAP BEHAVIOR 

We study in this section the reduction of of from T*^^ due to the fluctuations of the collective modes, using 
the Ginzburg-Landau free energy derived earlier. As a warm-up, as well as for the purpose of later comparison, 
we study the same effect in a weak coupling BCS superconductor with shprt range ((5-function) interaction, whose 
Ginzburg-Landau free energy is the familiar 0{N) 0* theory with N = 2-!E3 

Se = l3 [ dr(~-^\\/^/j\^ + a\ij\^ + b\ijA , (5.1) 



where in 3D we have 



4me 



/7C(3)meA:i?ei^y^^ 
^('•^ = [ l2.^{TMFy ) ^W' (5-2) 



7C(3)eF 
9^4(rMf)2 

14C(3)mefcFe|^ 



(5.4) 



and rUe is the effective mass of the electron. 

In the mean field theory, one neglects the quartic term in (^^), and Tc is determined by the point where a turns 
negative, which is nothing but T*^^. The fluctuation effects oi_jbhe quartic term may be studied using the self- 
consistent field approximation that is exact in the limit N oorB Within this approximation, one replaces IV'I'* by 
4(|V'P)|V'P; and requires that the following self-consistency equation is satisfied: 

d^l') = 17 T^^M-^.) = 7^3/ ^'k^^— — . (5.5) 

At T = we ought to have a -t- 46(|7/'p) = 0. Thus the equation that determines T^ becomes 



AmckBTc f d^k 



^b^^ I — = 0, (5.6) 



hujo /vp is an appropriate high momentum (ultraviolet) cutoff for k. From this we obtain 



- Tc knTchuj 



« 1 (5.7) 



J^MF 

for weak coupling superconductors, whence it is clear that the mean field value for Tc is extremely accurate. 

In the following we demonstrate that the situation is very different, and Tc becomes significantly lower than T*^^ 
when L ^ ^, due to thermal fluctuations of the low energy collectiv e mo des discussed in previous sections. We 
use the Ginzburg-Landau theory appropriate for this situation, Eq. ( 4.29| ), and use the same self-consistent field 



approximation as above by replacing |A(R, k)|^ with 4(| A(R, k)p) | A(R, k)p, after which the reduced action takes a 
quadratic form: 
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[A, = ^ ^ y dR I - CiP, k) + ?7(fc)]| A(R, k)|2 + Bal VkA(R, k)p + • VA(R, k) 



(5.8) 



where 

;7(fc)=4;7(fc)(|A(R,k)|2). (5.9) 

We would like to formally diagonalise this quadratic form. To do that, we take advantage of the translation invariance 
of the action, and express the reduced action in terms of the Fourier transform of A(R, k), A(Q, k): 



k,Q 



|A(Q,k)|2+B2|VkA(Q,k)|n . (5.10) 



2mk 

Note that modes with different Q's decouple. To proceed further, we introduce eigen modes ipmn{Q, k) which satisfy 



Bi - C(/3, k) + U{k) + !^-S)l _ B2VI 

2nik 



V'mn(Q,k) = i;„„((5)-0mn(Q,k), (5.11) 



and the normalization condition: 

k 

Here m and n are quantum numbers to be specified later. Expanding A(Q,k) in terms of ?/;mn(Q,k): 

A(Q, k) = ^ a^niQHmniQ, k), (5.13) 

mn 

we bring the reduced action to diagonal form: 

Sr[A,A] = ^Y1 Emn{Q)\aran{Q)\^- (5.14) 



A 

Qmn 



The self-consistent equation now becomes 



4t/(fc) V- , , .r^ , m2/, ,^m2x _ 4{/(fc) |Vw(Q,k)| 

A' 

At T = Tc, we have the lowest eigenvalue 



C/(fc) =4L/(fc)(|A(R,k)|^) = ^ E |V.™«(Q,k)|Ma_(Q)|^) = ^ E eI^Q) ' ^^''^^ 

Qmn Qmn 



EoaiQ = 0) = 0; (5.16) 

This poses a self-consistent conditi on th at determines Tc. 

The Schroedinger-like equation (5.11) is quite difficult to solve in general, one of the reason being the isotropy in k 
space is lost due to the presence of the term proportional to (k • Q)^. In order to proceed, we need to make a number 
of further approximations and assumptions. Firstly, we average (k • Q)^ along different directions of k and replace it 
by (5^/2, so that isotropy in k space is restored and all modes may be labeled by an "angular momentum" quantum 
number m. We also notice that —C{P,k) has its minimum at fc = kp^ and goes to zero rapidly for > 1//3. Hence 
we may approximate it by a "square" well, with depth — C(/3, fc^) = f and width i^jjj^- U{k) (and therefore U{k)) 

has a similar structure, and we can also approximate it by a step-like structure, with height Ukp = and the same 
width . Anticipating that the "low-energy" (E) modes will be localized in the step well, we may approximate 
by a constant . With these simplifications, Eq. ( 5.11 ) reduces to that of the Schroedinger equation of a particle 
confined to a step potential well. Imposing periodic boundary condition at the ends of the welljs^ we obtain 



V'mn(Q,k) = /^e*[2™e,+„.(/c-/c,)/../3]^ ^5^^^) 
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and at T = Tc (using ( |5.16D ) 



"64 ' ^^^P 



^mn(Q) = + B2(tX + TT^ri^^r )• (5.18) 



Combining Eqs. ( 3.15| ) and (5.18), we obtain the equation that determines Tc'. 



We notice that for the case of m = and n = 0, the integral is logarithmically divergent at both infrared and 
ultraviolet. As discussed earlier, there is always an ultraviolet cutoff A„.„ « Ec/vp- The infrared divergence is a 
signature of the fact that 2D is the lower dimensionality for ordering of this model (we do not go into details of 
the Kosterlitz-Thouless picture for a true phase transition for N = 2 here). An infrared cutoff A^^ is provided by 
invoking the quasi-2D nature of all real systems, which eventually crossover to 3D at sufficiently long length scales. 
Assuming A^^ is sufficiently large compared to the energies (measured in proper units) of the modes whose fluctuations 
contribute significantly to the reduction of Tc, we obtain 

^ 4M ^ A„„ 
3vpkp Air 

where M is the number of modes contributing in the sum of m and n. 

To determine M, we note that the self-consistent potential (for Q — 0; finite Q only adds a constant to it) is zero 
inside the well, and Bi outside it. Thus in order for the previous approximations for mode solutions and "energies" 
to be valid, we need to have Emn{Q = 0) < Bi. Summing up the number of these modes, we obtain 



M- kF\/Bl/B2 = kpL. (5.21) 

We thus find 



VP T^f^ 

It clearly goes to zero as L — > cx).0 We would like to emphasize that here we have made a number of crude ap- 
proximations in the calculation of Tc, thus the dependence of Tc on the range L may not be quantitatively reliable. 
However it is quite clear that the approximations we made tend to underestimate the importance of the fluctuation 
effects; thus the qualitative conclusion that Tc — > as i — > oo must hold. 



VI. SUMMARY AND DISCUSSION 



In the bulk of the paper, we have taken the attitude that the model we study here, namely a system with long 
(compared to the coherence length) but finite range pairing interaction, is a theoretical model that is interesting in 
its own right, and worked out some unusual properties of this model, with emphasis on those properties that are 
qualitatively different from those of the standard weak coupling BCS superconductors stabilized by a short-range 
pairing potential. Our most interesting finding is that in this model the transition temperature Tc is controlled by 
thermal fluctuations of collective modes, and can be significantly lower than the quasiparticle gap A or the mean- 
field transition temperature T*^^; as a consequence in the temperature range Tc < T < Tc '^ the systerp-jexhibits 
pseudogap behavior as the electrons are still paired while there is no superconducting long-range order.ES In this 
section we attempt to make contact between our results and the phenomenology of cuprate superconductors, discuss 
the relation between our model and existing theoretical work on the pseudogap behavior, point out the limitations of 
our model as well as of our analysis, and indicate some natural extensions and directions for future study. 

One of the motivations of the present study is the observation that the coherence length ^ is much shorter in the 
cuprates than in conventional superconductors. We are, however, by no means the first to suggest that a short ^ can 
lead to behavior qualitatively different from weak coueling BCS theory. Since tha pady days of high TcirffiUfwing 
the work of Leggett and Nozieres and Schmitt-Rink,E3 Randeria and coworkers, &E3eI as well as others,c2lEil have 
argued that the short coherence length may bring the cuprates to a regime that is intermediate between the weak- 
coupling BCS limit and the Bose-Einstein Condensation (BEC) limit of Cooper pairs. The latter case is realized if 
^ is much shorter than the inter-particle spacing so that kp£, ^ 1; in terms of energy scales, that corresponds to 
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the case A ^ Ep. In this case the Cooper pairs are so closely bound that they hardly overlap, and may be viewed 
as point-like hard-core bosons at low energies, while the transition temperature Tc is essentially their Bose-Einstein 
condensation temperature which is much lower than and unrelated to the pairing energy A. For Tc < T < A, the 
electrons remain paired yet there is no long-range superfluid order, hence pseudogap behavior. In terms of low-energy 
excitations responsible for destroying superconductivity, in the BEC limit it is the linear Goldstonc mode (assuming 
no Coulomb interaction), while the quasiparticle excitations (broken pairs) cost too much energy to have any effect 
on Tc- Put differently, the thermal fluctuations of this Goldstonc mode are the classical phase fluctuations of the. 
superconducting order parameter that control the transition in this limit. (We should note that Emery and KivelsonEj 
have argued that classical phase fluctuations are the physics of the pseudogap regime in the underdoped materials, 
but their point of departure is logically distinct, deriving from a small zero temperature superfluid stiffness that 
is connected to the physics of a doped Mott insulator.) In the cuprates, kp£, ^ 10, and A is still a small energy 
compared to Ep (although the difference is not nearly as overwhelming as in conventional superconductors); we are 
thus still somewhat distant from the BEC limit. The interesting new feature of the model studied here is that by 
having ^ much smaller than the range of pairing interaction L, we can get the pseudogap behavior while staying the 
weak coupling regime (in the sense kp£^ ^ 1 and A/Ep ^1). In this case the Goldstonc mode is unable to drive Tc 
below the scale of A by itself; it needs all the help from the other collective modes supported by this model. The 
low-energy spectra of the weak-coupling BCS superconductor, the BEC superconductor, and the present model are 
summarized schematically in Fig. l|. While there are clearly qualitative differences between the BEC (as well as 
phase fluctuation) picture and the present model, they share the common spirit that T^, is determined by collective 
modes instead of quasiparticle excitations. It is also worth mentioning that in real systems with long-range Coulomb 
interaction, the energy of the Goldstonc mode is expected to be pushed up to the plasmon frequency; the importance 
of this fact to the thermal fluctuations of this mode (or classical phase fluctuations) is still under discussion.Ej On 
the other hand the Coulomb interaction is not expected to significantly affect the spectra of the (gapped) exciton-like 
modes discussed here. This is due to the different nature of the Goldstonc mode and exciton modes: the former 
is a consequence of the broken gauge symmetry, and introducing the Coulomb interaction (or, more generally, the 
electromagnetic interaction) converts the Goldstonc mechanism of broken symmetry to the Higgs mechanism. The 
exciton modes, on the other hand, are the bound states formed by quasiparticle pairs due to the residual attractive 
interaction between quasiparticles not included in the mean field approximation. We therefore believe the presence 
of Coulomb interaction does not affect the fluctuation physics discussed here significantly. 

As pointed out already in the Introduction, in some theoretical models for cuprate superconductivity, the-cange 
of the interaction that gives rise to Cooper pairing can be very long. In the interlayer pair hopping model,allJ it is 
assumed that pairing is induced by a pair hopping term in the Hamiltonian, that is diagonal in momentum space. 
Fourier-transforming to real space, this corresponds to an infinite range hopping term for Cooper pairs, corresponding 
(loosely) to the L ^ oo limit of the model we study here (the fact that the hopping terar is off-diagonal in layer 
index is of no qualitative consequence). In that specific form, it has already been showncB that the model may be 
solved exactly in the absence of any other in-plane paring interaction, the transition temperature is zero, and there 
is pseudogap behavior at low temperatures. Our results are in agreement with this observation. The model we use 
here, however, is more general and versatile, and in particular enables one to address how the pseudogap behavior 
develops as the range L increases, and how Tc approaches zero as L — s- oo. l-l 

It is equally interesting to scrutinize the results obtained here in the context of the spin fluctuation theoryu ofxiiprate 
superconductivity. In this model the range of the interaction L is essentially the spin-spin correlation lengthHj Thus 
L is of order lattice spacing in the overdoped region of the phase diagram, increases as the doping level x decreases, 
becomes much longer than the lattice spacing in the underdoped region, and eventually diverges upon approaching 
the antiferromagnetic phase boundary at very low doping. The pseudogap behavior is observed in the underdoped 
regime, where Tc — s- while the gap (the maximum of the c?-wave gap measured at very low T by, say, photoemission) 
slowly increases as x decreases. This behavior is certainly consistent with our findings here; the reduction of Tc and 
its departure from the gap is due to the increase of L, while the size of gap, which in our model saturates at the depth 
of the potential well Vo for large L, would be set by the scale of the near neighbor spin coupling strength J. Hence 
at a very crude level, we can qualitatively account for the phenomenology of the underdoped cuprates by combining 
our results with the spin-fluctuation model. 

On the experimental side, some circumstantial evidence for both the importance of the longer range part (beyond 
near-neighbor)-pf the pairing interaction, and the possible relevance of our results, exists, i) Recent photoemission 
measurements^ of gap anisotropy have found deviations from the standard 

Ak C!C cos fca; — cos fcy (6-1) 

dependence of the d,j.2_y2 order parameter in tka underdoped region, with the deviation increasing with the decreas- 
ing doping level. This deviation is interpretedl23 as due to longer-ranged pairing interaction, as a nearest-neighbor 



16 



attraction leads to Eq. (6.1). We consider this to be (somewhat indirect) evidence that the range of the-pairing in- 
teraction increases with decreasing doping level in underdoped cuprates. (ii) It has been recently noticedE3 that there 
is a strong correlation between Tc and JPfak width of the normal state spin susceptibility x((7,a'), in some cuprates. 
Specifically, Balatsky and and BourgesEj found that in YBCOi27i Etnd La2u compounds, Tc oc 5q, where 5q is the 
width of ImxiQ,'^) ^^"^ ^- This lead the authors to conclude that "antifcrromagnetism is likely responsible for the 
high Tc superconducting mechanism". If so, one is lead to the conclusion that Tc oc I /L, where L ^ 1/Sq is the range 
of th e pai ring interaction mediated by the normal state antiferromagnetic spin fluctuations, in agreement with our 
Eq. ( 5.22 )! However, given the facts that our model is greatly oversimplified as far as the cuprates go (see below) and 
our treatment of the fluctuation effects is still preliminary (also below), we consider this agreement to be fortuitous 
at this point. Nevertheless, the findings of Ref. |3^ do indicate the importance of thje-xange of interaction on Tc, and 
are in qualitative agreement with our results, (iii) Experimentally, it has been foundcj that upon cooling underdoped 
cuprate samples, the temperature T* at which the pseudogap opens up depends on the location in momentum space; 
T* is highest near the superconducting gap maxima, while lowest near the gap nodes; it thus suggests that T* is 
a "local" property in momentum space. This can be understood quite naturally by invoking a long range pairing 



interaction as studied here, since such an interaction is localized in momentum space. As we have shown in section II] , 
one can define a momentum-dependent mean field transition temperature, ^^^^^(ek), which increases monotonically 
with Ak; it also sets the temperature T* ~ Tc'^{e-k) below which a local gap is opened up and the pseudogap sets in. 

We must emphasize that the contact between our model and the cuprate physics made above is tentative, and in 
its present form this model can only be viewed as an interesting toy model that gives rise to pseudogap behavior. 
In the following we discuss some limitations of the model as well as our treatment, and some natural directions for 
extensions. 

The ground state of our model is a fully-gapped s-wave superconductor, while the cuprates (at least most of them) 
are known to have a d-wave order parameter. The most important, difference between the two is that the latter 
supports gapless nodal quasiparticles. It has already been suggestedlj that the thermally excited quasiparticles may 
be responsible for emptying the superfluid stiffness, setting the scale of Tc, and giving rise to pseudogap behavior 
in underdoped cuprates. This important piece of physics is missing in our model. On the other hand, as long as 
collective modes of the order parameter are concerned, our method can be generalized to d-wave (or other unconven- 
tional superconductors) fairly easily. The key feature of a long (but finite) range pairing interaction is its locality in 
momentum space, i.e., Vk is sharply peaked in k space. This allows a gradient expansion in momentum space for the 
order parameter. In an s-wave superconductor, Vk is sharply peaked at k = 0, while for a d-wave superconductor 
appropriate for cuprates, one would choose a Vk that is sharply peaked at a wave vector Q which is at or near (tt, tt). 
This, however, is not going to make any qualitative difference to the analysis made in this paper. At a more detailed 
level, two cases need to be distinguished from each other, (i) There is Fermi surface nesting and Q is exactly or 
very close to the nesting wave vector. In this case the entire Fermi surface participates in pairing actively and the 
analysis carried out here, based on the gradient expansion of the superconducting order parameter both along and 
perpendicular to the Fermi surface, carries through straightforwardly, (ii) There is no Fermi surface nesting, or Q is 
not close to the nesting wave vector. In this case only certain "hot spots" at the Fermi surfaces that are connected 
by Q participate in the pairing actively. In this case one can still develop a gradient expansion within the hot spots. 
In either case the qualitative [features of our results are expected to be robust. We note that in a recent study of the 
spin-fermion-hot spot model,E3 it has been conjectured that in the limit of infinite spin-spin correlation length, the 
transition temperature obtained from solving the Eliashberg equation is only the onset of pseudogap behavior, while 
the real Tc vanishes in that limit. This conjecture is in agreement with our results, as the spin-spin correlation length 
corresponds to the range of interaction L in our model. In the present paper, we have developed a systematic way to 
study the physics of pseudogap behavior due to large L, and shown explicitly that Tc ^ as L — > oo. 

Our analysis of the thermal fluctuations of the collective modes is based on a Ginzburg-Landau free energy func- 
tional, which we derive using a functional integration formalism and an expansion in power series of the magnitude 
of the order parameter. Strictly speaking this power series expansion is valid only in the vicinity of T*^^, where the 
amplitude of the order parameter just starts to develop; aX T = Tc T^"^^ , the amplitude of the order parameter is 
big and such an expansion is no longer appropriate. What is also missing are the contributions from the fluctuations 
of components of the order parameter with finite Matsubara frequencies; we neglected them on the ground that we are 
primary interested effects of thermal fluctuations; however at low T these components are important to the physics, 
especially if one is also interested in the quantum fluctuations of the order parameter. Even with these simplifications, 
the resultant free energy functional is quite complicated, and we need to introduce approximations in the one-loop 
calculation of Tc, like neglecting the anisotropy in the internal (relative) space of the order parameter when the center 
of mass carries a finite momentum. It is quite possible that one can do a better job in analyzing the fluctuation, 
especially in the limit of L/^ oo. It is worth noting, however, that all our approximations tend to underestimate 
the effects of fluctuations, and thus our basic conclusion that Tc <^ Tc^^ ~ A when L/^ ^ 1, is clearly vahd. 

Despite the various disclaimers made above, we hope the model and the crossover introduced here will serve as a 
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different paradigm of a non-BCS transition to a superconducting state, which can exhibit pseudogap behavior even 
without leaving the weak coupHng regime. Its relevance to cuprate physics is not completely clear at present, but 
there are encouraging signs that it is worth further pursuit. 
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FIG. 1. Schematic illustration of the neutral excitation spectra of three diilerent types of superconductors, in the absence 

of Coulomb interactions. The region shaded by dashed lines with a gap 2A is the quasiparticle pair continuum, common to 
all three cases; the solid lines stand for collective modes. For a weak-coupling BCS superconductor with short-range pairing 
interaction (part (a): i; ^ l/kp and ^ S> L, where ^ is the coherence length, kp is the Fermi wave vector, and L is the range 
of the pairing interaction), the only collective excitation is the linear Goldstone mode, whose velocity {v ~ vf) is very big as 
compared to the scale of the gap: vkr ^ 2A, thus contributes little to the thermodynamics and hence the determination of 
Tc. In the BEG case (part (c); ^ ^ l/kp), again the linear Goldstone mode is the only low-energy collective mode, but due to 
the reversed energy scale; vkp <^ 2A, it dominates the thermodynamics and sets the scale for Tc- The situation studied in this 
paper, as illustrated in part (b), is in some sense between these two extremes, namely we have ^ ^ l/kp as in weak-coupling 
BCS, but in the meantime ^ <^ L. Here the Goldstone mode is very steep and hence contributes little to the thermodynamics; 
however the the combined effect of all the low-energy modes, stabilized by the condition ^ -C i, dominates the thermodynamics 
and determines Tc. 
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